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Abstract. A length- TV spin chain with the y/N(— v)-th neighbor interaction is 
identical to a two-dimensional (d = 2) model under the screw-boundary (SB) condition. 
The SB condition provides a flexible scheme to construct a d > 2 cluster from 
an arbitrary number of spins; the numerical diagonalization combined with the SB 
condition admits a potential applicability to a class of systems intractable with the 
quantum Monte Carlo method due to the negative-sign problem. However, the 
simulation results suffer from characteristic finite-size corrections inherent in SB. In 
order to suppress these corrections, we adjust the screw pitch v(N) so as to minimize 
the excitation gap for each N. This idea is adapted to the transverse-field Ising model 
on the triangular lattice with N < 32 spins. As a demonstration, the correlation-length 
critical exponent v is analyzed rather in detail. 
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1. Introduction 

A length- N spin chain with the y/N (= t>)-th neighbor interaction is identical to a two- 
dimensional (d = 2) model under the screw-boundary condition; a schematic drawing 
is presented in Fig. [TJ The screw-boundary condition provides a flexible scheme to 
construct a d > 2 cluster from an arbitrary number of spins N. With the aid of Novotny's 
method [TJ [2] , one is able to implement the screw-boundary condition systematically; 
details are explicated afterward. Because the quantum Monte Carlo method is 
inapplicable to a class of systems such as the frustrated quantum magnets (negative- 
sign problem), the numerical diagonalization combined with Novotny's method might 
provide valuable information as to such an intriguing subject. 

As a matter of fact, the screw-boundary condition was adapted to the S = 1 
triangular magnet for N < 20 [3]; tuning the spatial anisotropy and the biquadratic 
interaction, we are able to realize the deconfined criticality jl] between the resonating- 
valence-bond (RVB) and magnetic phases. The present scheme (modification) is readily 
applicable to this problem. Moreover, a recent reexamination of the transverse-field 
Ising model in d — 2 and 3 via the real-space renormalization group revealed that a 
procedure, the so-called "exact invariance under renormalization," yields a substantially 
improved critical exponent even for an analytically tractable system size; a relevance to 
the present treatment is addressed in the last section. 

In this paper, we make an attempt to reduce the finite-size corrections inherent in 
the screw-boundary condition. As a matter of fact, the simulation results exhibit bumpy 
finite-size deviations depending on the condition whether the pitch v{= y/N) is close to 
an integral number or not; actually, in Ref. [Tj, such a singular case is treated separately. 
The screw pitch v admits a continuous variation, and there is no reason to stick to a 
fixed pitch v = y/N. We adjust y/N) so as to minimize the excitation gap AE 
for each system size N. As a demonstration, we adapt this idea to the (ferromagnetic) 
transverse-field Ising model on the triangular lattice for iV < 32. In order to demonstrate 
the performance of this approach, we analyze the criticality (correlation-length critical 
exponent v) rather in detail; the simulation results are compared with those of the fixed 
screw pitch v = y/N. 

The rest of this paper is organized as follows. In Section [2J we explicate the 
simulation algorithm, placing an emphasis on the scheme how we fix the screw pitch to 
an optimal value. In Section El we show the numerical results. As a demonstration, we 
make an analysis of the critical exponent v. In Section HJ we present the summary and 
discussions. 

2. Simulation algorithm for the two-dimensional transverse-field Ising 
model under the screw-boundary condition: Novotny's method 

We employ Novotny's method in order to implement the screw-boundary condition 
[TJ [2] . For the sake of selfconsistency, in this section, we present an explicit simulation 
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algorithm. 

Before commencing an explanation of the technical details, we sketch a basic idea 
of Novotny's method; see Fig. [TJ We implement the screw-boundary condition for a 
finite cluster with N spins. We place an S = 1/2 spin (Pauli matrix Bj) at each lattice 
point i(< N). Basically, the spins constitute a one-dimensional (d = 1) structure. The 
dimensionality is lifted to d = 2 by the long-range interactions over the y/~N(= v)-th- 
neighbor distances; owing to the long-range interaction, the N spins form a y/~N x \/N 
rectangular network effectively. (The significant point is that the screw pitch v ~ \/N 
is not necessarily an integral number, and can even afford to vary continuously.) 

According to Novotny [H [2], the long-range interactions are introduced 
systematically by the use of the translation operator P; see Eq. The operator 

P satisfies the formula 

P\ai, a 2 ,...,a N ) = \a N ,a 1 , . . . , erjv-i)- (1) 
Here, the Hilbert-space bases o~ 2 , ■ ■ ■ , o"v)} (&i = ±1) diagonalize the operator of; 

^IW> = ^ilW>- (2) 

2.1. Details of Novotny's method 

We formulate the above idea explicitly. We adapt Novotny's method to the transverse- 
field Ising model on the triangular lattice. The Hamiltonian with a screw pitch u(« y/N) 
is given by 

N 

-rj>*, (3) 

i=l 

with the transverse magnetic field V. Hereafter, we consider the exchange interaction J 
as a unit of energy (J = 1). The i>-th neighbor interaction H(v) is a diagonal matrix, 
whose diagonal element is given by 

H WlU ^) = (WI^(«)IW> = (fo}|TF»|fo}>. (4) 

The insertion of P v is a key ingredient to introduce the i>-th neighbor interaction. Here, 
the matrix T denotes the exchange interaction between {o"j} and {t{\\ 

N 

({a t }\T\{T. l }) = Y,°kT k . (5) 

k=l 

Afterward, we search for an optimal value of the screw pitch u(ps y/N). 

The above formulae complete the formal basis of our simulation scheme. We 
diagonalize the Hamiltonian matrix ()3]) for N < 32 spins numerically. In the practical 
numerical calculation, however, a number of formulas may be of use; see the Appendices 
ofRefs. jSHZ]. 
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2.2. Search for an optimal screw pitch v 

An optimal value of the screw pitch v is determined as follows. In Fig. |2j we plot the 
(first) excitation gap AE for the parameter a = v 2 /N, and the system size, N = (+) 
28 and (x) 32. Here, we fix the transverse magnetic field to F — 4.8, which is close 
to the critical point (separating the paramagnetic and ferromagnetic phases); see Fig. 
[3j (An underlying idea behind this analysis is presented afterward.) From Fig. [21 we 
see that the excitation gap takes a minimum value beside a = 1; note that so far, the 
screw-boundary-condition (Novotny method) simulation has been performed at a = 1 
(f = \/N). We adjust v(N) so as to minimize the excitation gap for each N. As a 
result, we arrive at an optimal screw pitch (v = V Na) 

a = 1.12 , 0.98 , 0.81 , 0.68 , 1.30 , and 1.06 , (6) 

for N = 12, 16, ... , 32, respectively. (The series of results may indicate a regularity such 
that v^V-dependence of a is like sawtooth.) The optimal screw pitch differs from the 
conventional one [Tj, |2] 

v = VN. (7) 

A comparison between Eqs. (jSJ) and (jTJ) will be made in order to elucidate the 
performance of the present approach. 

The underlying idea behind the above screw-pitch optimization (jBJ) is as follows. In 
general, the excitation gap is proportional to reciprocal correlation length £ oc 1/AE. 
On the one hand, according to the scaling hypothesis, namely, at the critical point, the 
correlation length develops up to the linear dimension of the cluster L ~ £; in other 
words, the correlation length is bounded by the system size, and it should extend fully 
for the equilateral cluster. Hence, as the energy gap gets minimized, the (effective) shape 
of the cluster restores its spatial isotropy, rendering the finite-size behavior smoother. 

3. Numerical results 

Setting the screw pitch to an optimal value (jHJ), we turn to the analysis of the exponent 
v rather in detail. In Fig. [31 we plot the scaled excitation gap for various T and 
AT = 12, 16, ... , 32; the system size L is given by L = y/N, because N spins constitute a 
rectangular cluster. According to the scaling hypothesis, the intersection point indicates, 
a location of the critical point; for details of the scaling theory, we refer readers to a 
comprehensive review article jH]- We see that a phase transition occurs at V « 4.8. 

In Fig. HI we plot the approximate critical point T c (Li,L 2 ) for [2/(Li + L 2 )} 3 with 
12 < Ni < N 2 < 32 (Li j2 = y/Nip). [The abscissa scale 1/L 3 is based on the scaling 
theory, which states that in principle, the corrections of r c should behave like 1/ L l l u+U> 
with 1/u + u = 1.5868(3) + 0.821(5) j9]. There may also exist subdominant corrections, 
and the abscissa scale is set to 1/L 3 .} Here, the symbol (+) denotes the simulation 
result with the optimal screw pitch adjusted to Eq. (jSJ), whereas the data of (□) were 
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calculated with the screw pitch fixed to Eq. fl7|) tentatively. The approximate critical 
point satisfies the fixed-point relation with respect to LAE 

Li AE (Li ) | r=r c (ii ,L a ) = L 2 AE{L 2 )\ T=rc{LlM) , (8) 

for a pair of system sizes (L 1 ,L 2 ). The data with the conventional scheme [symbol 
(□)] exhibit a pronounced bumpy finite-size behavior, at least, as compared to those 
of (+). The bottom of the (□) data corresponds to iV « 5 2 (1/5 3 = 0.008); for such 
a quadratic system size, the data suffer from an irregularity pQ intrinsic to the screw- 
boundary condition. Such irregularities appear to be remedied by optimizing v [symbol 
(+)] satisfactorily. The extrapolated critical point rd^-^ is no longer utilized in the 
subsequent analyses; rather, the finite-size critical point r c (Li,L 2 ) is used. Hence, we 
do not go into further details as to the reliability of T c . 

In Fig. [5j we present the approximate correlation-length critical exponent 

lAT M = HLi/L 2 ) 

1 u 2) ln{9 r [ J L 1 AE(L 1 )]/5 r [L 2 A J E;(L 2 )]}| r=rc(Xl)X2) ' 1 ) 

for [2/(Li + L 2 )] 2 with 12 < N x < N 2 < 32 (L h2 = y/N^). [The abscissa scale 1/L 2 
is based on the scaling theory, which states that in principle, the corrections of critical 
indices should behave like 1/L W with u = 0.821(5) [9J. Afterward, we make a slight 
variation of the abscissa scale so as to appreciate possible subdominant contributions 
properly] The symbols, (+) and (□), denote the results with the optimal [Eq. (jB])] 
and fixed [Eq. ((7j)] screw pitches, respectively. Again, the former result (+) appears to 
exhibit reduced finite-size corrections as compared to the latter one (□). 

We are able to calculate the critical exponent, based on the Binder parameter 

u - l ~lW?' < 10 > 

with the magnetization M = Yli=i a t an d the ground-state expectation (...). In Fig. 
[5j we present the approximate critical exponent 

v{ T M = HLi/L 2 ) 

1 u 2) HWiLj/drUilvWr^.frM)' 
with 12 < Ni < N 2 < 32 and the approximate critical point r c (Li,L 2 ) satisfying the 
fixed-point relation 

C/(Li)|r=r c (i,i,z,2) = ^(^2)|r=r c (i,i,i 2 )- (12) 

The symbols, (x) and (o), denote the results calculated with the optimal [Eq. (JSJ)] 
and fixed [Eq. ((Tj)] screw pitches, respectively. Again, the former result (x) appears to 
exhibit reduced finite-size deviations. 

Apart from an improvement due to the optimal screw pitch, in Fig. [5j it would 
be noteworthy that the data of the symbols, (+) and (x), are complementary. That 
is, these results, obtained independently via the scaled energy gap [Eq. dHJ)] and the 
Binder parameter [Eq. (fTTT)]. provide the upper and lower bounds for u, respectively. 
Encouraged by this observation, we carry out an extrapolation to the thermodynamic 
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limit for v. An application of the least-squares fit to the data (+) yields an estimate 
v = 0.635(4) in the thermodynamic limit L — > oo. Similarly, the result (x) yields 
v = 0.619(5). An appreciable systematic deviation seems to exist. Replacing the 
abscissa scale of Fig. [5]with [2/(L 1 + L 2 )] 1 ' 8 , we arrive at the estimates v = 0.628(4) and 
v = 0.624(6) for the data (+) and (x), respectively; the deviation is now within the error 
margins. These results may lead to an estimate v = 0.626(10) with the error margin 
covering both systematic and insystematic (statistical) errors. As a comparison, we refer 
to a recent Monte Carlo simulation result v = 0.63020(12) [9] for the three-dimensional 
(classical) Ising model. This result is consistent with ours within the error margins, 
suggesting that possible systematic errors are appreciated properly. (As mentioned 
below, an application of the perfect-action technique might admit a suppression of the 
error margin.) 

Last, we address a remark. In order to demonstrate a reliability of the present 
approach, a comparison with an existing result would be of use; however, the Ising 
model on the triangular lattice has not been studied very thoroughly. In this respect, 
we refer readers to Ref. [10] (see Fig. 4), where the ordinary Ising model was simulated 
under the screw-boundary condition with v = y/N (a = 1); the critical point appears to 
be close to the existing result. An aim of this paper is to reduce the irregular finite-size 
corrections of Fig. 4 in Ref. [10] . 

4. Summary and discussions 

An attempt to suppress the finite-size corrections inherent in the screw-boundary 
condition (Fig. [1]) was made for the transverse-field Ising model on the triangular 
lattice. The screw pitch can afford to vary continuously, and there is no reason to stick 
to a fixed pitch v = \/~N. Setting the transverse magnetic field to a critical-point value 
T = 4.8, we scrutinize v-dependence of AE (Fig. [2]). Minimizing AE, we determine an 
optimal value of v(N) for each N, Eq. ([6]), where an effective shape of the cluster should 
restore its spatial isotropy. Actually, the numerical results with the optimal screw pitch 
for iV < 32 exhibit suppressed finite-size corrections, rendering the convergence to the 
thermodynamic limit smoother (Figs. H] and |5|; As a byproduct, the data obtained 
from the energy gap [Eq. (Q] and the Binder parameter [Eq. (fTTj) ] provide the upper 
and lower bounds for v, respectively. Encouraged by this observation, we make an 
extrapolation to L — > oo to obtain an estimate v = 0.626(10) via replacing the abscissa 
scale with 1/L 18 . 

We mention a number of remarks. First, as mentioned in Introduction, the 
transverse-field Ising model in d = 2 and 3 was reexamined with the analytical real-space 
renormalization group [5|; even for such an analytical tractable system size, detailed 
scrutiny of the scaling procedure affords us valuable information as to the criticality. 
To be specific, one may eliminate the residual finite-size corrections (~ 1/L 18 ) by 
extending the interaction parameters and tuning them through the perfect-action 
technique [TU [121 E3 ELH HSJ US] ; this technique has been developed in the context of the 
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Figure 1. A schematic drawing of the spin cluster for the transverse-field Ising model 
Q is presented. We implement the screw-boundary condition with the screw pitch v, 
based on Novotny's method. As indicated above, the spins constitute a one-dimensional 
(d = 1) alignment {di} (i = 1, 2, . . . , N), and the dimensionality is lifted to d — 2 by the 
w(ss v / A)-th- neighbor interactions. Full details of the simulation scheme are presented 
in the text. 

lattice field theory, for which the continuum (thermodynamic) limit has to be undertaken 
properly and efficiently Second, the present scheme (modification) is readily applicable 
to the quantum two-dimensional magnet [3], which exhibits a deconfined criticality 
(phase transition between the RVB and magnetic phases). Because of the negative-sign 
problem, the quantum Monte Carlo method is not quite efficient to realize the RVB 
phase, and the numerical diagonalization combined with the modified screw-boundary 
condition would provide a promising candidate so as to survey such an intriguing issue. 

References 

[1] M.A. Novotny, J. Appl. Phys. 67 (1990) 5448. 

[2] M.A. Novotny, Phys. Rev. B 46 (1992) 2939. 

[3] Y. Nishiyama, Phys. Rev. B 83 (2011) 054417. 

[4] R.R.P. Singh, Physics 3 (2010) 35. 

[5] R. Miyazaki, H. Nishimori, and G. Ortiz, Phys. Rev. E 83 (2011) 051103. 

[6] Y. Nishiyama, Phys. Rev. E 75 (2007) 051116. 

[7] Y. Nishiyama, Nucl. Phys. B 832 (2010) 605. 

[8] A. Pelissetto and E. Vicari, Phys. Rep. 368 (2002) 549. 

[9] Youjin Deng and H. W. J. Blote, Phys. Rev. E 68 (2003) 036125. 

[10] Y. Nishiyama, Phys. Rev. E 74 (2006) 016120. 

[11] J.H. Chen, M.E. Fisher, and B.C. Nickel, Phys. Rev. Lett. 48 (1982) 630. 

[12] K. Symanzik, Nucl. Phys. B 226 (1983) 187. 

[13] P. Hasenfratz and F. Niedermayer, Nucl. Phys. B 414 (1994) 785. 

[14] H.W.J. Blote, E. Luijten, and J.R. Heringa, J. Phys. A: Math. Gen. 28 (1995) 6289. 



An optimization of the screw pitch 



8 



LU 

< 




1.28 ' 1 1 1 1 1 1 

0.9 1 1.1 1.2 1.3 1.4 

a 

Figure 2. The excitation gap AE is presented for a parameter a = v 2 /N (v is the 
screw pitch), and the system size, N = (+) 28 and (x) 32; the transverse magnetic 
field is fixed to T — 4.8, which is close to a critical point (Fig. [3]). We determine an 
optimal screw pitch v(N) so as to minimize AE for each N. 




Figure 3. The scaled energy gap LAE is plotted for the transverse magnetic field T, 
and the system size N = (+) 12, (x) 16, (*) 20, (□) 24, (■) 28, and (o) 32 (L = y/N). 
Here, the screw pitch v is set to an optimal value [Eq. JBJ)]. We observe an onset of 
continuous phase transition around T m 4.8. 
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Figure 4. The approximate critical point r(Li,La) [Eq. (|SJ)] is plotted for 
[2/(Za + L 2 )} 3 with 12 < Ni < N 2 < 32 (L 1)2 = y/W[~2). The symbols (+) and 
(□) denote the data for the optimal [Eq. ([BJ] and fixed [Eq. ([7])] screw pitches, 
respectively. The data indicate that the optimal screw pitch suppresses the (bumpy) 
finite-size deviations. 
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Figure 5. The approximate correlation-length critical exponent v(Li,L 2 ) is plotted 
for [2/(Za + L 2 )} 2 with 12 < Nx < N 2 < 32 (i lj2 = ^/W^). The data of the 
symbols, (+) and (□), were calculated from the scaled energy gap [Eq. for the 
optimal [Eq. ([6])] and fixed [Eq. (0] screw pitches, respectively. Similarly, from the 
Binder parameter [Eq. (|11[) ]. the data of the symbols, (x) and (o), were calculated 
for the optimal [Eq. ((6j) ] and fixed [Eq. 0] screw pitches, respectively. The data of 
the optimized screw pitch, (+) and (x), exhibit suppressed finite-size deviations. An 
extrapolation to the thermodynamic limit is argued in the text. 
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